Wigner crystal in snaked nanochannels 
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We study properties of Wigner crystal in snaked nanochannels and show that they are charac- 
terized by conducting sliding phase at low charge densities and insulating pinned phase emerging 
above a certain critical charge density. The transition between these phases has a devil's staircase 
structure typical for the Aubry transition in dynamical maps and the Frenkel-Kontorova model. We 
discuss implications of this phenomenon for charge density waves in quasi-one-dimensional organic 
conductors and for supercapacitors in nanopore materials. 
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The Wigner crystal [1] appears in a great vari- 
ety of physical systems including electrons in two- 
dimensional semiconductor samples and one-dimensional 
(ID) nanowires (see review [2] and Refs. therein), elec- 
trons on a surface of liquid helium [3] , cold ions in radio- 
frequency traps [4] and dusty plasma in laboratory or in 
space [5] . Effects of Coulomb interactions are clearly seen 
experimentally in nanowires [6, 7] and carbon nanotubes 
[8]. Also interaction effects for electrons in microchan- 
nels on a surface of liquid helium have been recently 
observed experimentally [9]. In view of this remarkable 
progress it is interesting to investigate sliding and con- 
ducting properties of the Wigner crystal in wiggled or 
snaked nanochannels. The interest to such studies goes 
back to the Little suggestion [10] on electron conduction 
in long spine conjugated polymers where he proposed an 
approach for synthesizing of organic superconductors. A 
modern overview discussion of this important suggestion 
is given in [11]. A schematic image of electron transport 
in such organic molecules is shown in Fig. la. Accord- 
ing to this picture long molecules form wiggling channels 
which in principle can support electron transport along 
them. However, the Coulomb interactions between elec- 
trons are rather strong at such small scales and thus it 
is not obvious under what conditions a sliding of Wigner 
crystal along such channels is possible. This problem 
is related to the conduction properties of charge den- 
sity waves (CDW) (see e.g. reviews [11, 12]). To study 
this phenomenon we choose a simple model of ID snaked 
channel shown in Fig. Ib. There is no potential gradient 
along the channel but the channel walls are assumed to 
be very high so that electrons can move only along the 
channel. 

In addition, the properties of Wigner crystal in snaked 
nanochannels are also useful for understanding of mecha- 
nisms of charge storage in electrochemical capacitors, or 
supercapacitors, which start to have important industrial 
applications [13, 14]. In these systems, charged ions are 
stored in nanopores at the surface of the carbon-activated 
material which has enormously large capacitance C go- 
ing beyond the meanfield values given by the Helmholtz 



FIG. 1: (Color online) (a)A schematic image of the Little 
suggestion for electron transport in organic molecules (af- 
ter [10, 11]). (b) A schematic image of electron Wigner 
crystal with charges ei (points) sliding in a snaked sinu- 
soidal nanochannel, dashed lines show force directions be- 
tween nearby electrons. 



theory [15, 16]. At nanoscale the wiggling of pores is 
definitely present that makes our studies very timely. 

Due to sinusoidal channel wiggling the Wigner crystal 
moves in a certain effective periodic potential. The case 
of sliding of ID Wigner crystal in a periodic energy po- 
tential was analyzed in [17] having in mind an examplc of 
ion chains in optical lattices. It was shown there that this 
problem can be locally reduced to the Frenkel-Kontorova 
model for a particle spring chain in a periodic poten- 
tial [18] with particle positions described by the Chirikov 
standard map [19]. For a small amplitude of periodic po- 
tential the Wigner crystal with an incommensurate ion 
density can slide in an optical lattice but above a certain 
critical amplitude of potential the crystal is pinned by the 
lattice due to the Aubry analyticity breaking transition 
[20] . In the pinned phase the phonon Spectrum has a gap 
for long wave excitations so that this regime corresponds 
to an insulating phase. This situation corresponds to a 
dynamical spin glass phase with exponentially many sta- 
ble classical configurations being exponentially close to a 
ground state at a fixed electron density [17]. The Frenkel- 
Kontorova model is characterized by similar classical and 
quantum properties [21, 22]. At sufficiently large values 
of effective Planck constant a quantum instanton tunnel- 
ing between these quasidcgenerate configurations leads 
to a zero-temperature quantum phase transition due to 
which the quantum Wigner crystal becomes conducting 
[17]. In the following we show that the main elements of 
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FIG. 2: (Color online) Hull function s = h(x) (a,b) and 
phonon Spectrum u(k/N) (c,d) for incommensurate electron 
densities v = N/L = 239/233 (a,c) and v = N/L = 244/233 
(b,d). Here a = 1.2. 



this physical picture remain valid for the Wigner crystal 
in snaked nanochannels which are however characterized 
by enormously sharp changes of conducting properties. 

We start our analysis from the case of classical elec- 
trons with Coulomb interactions moving in a snaked 
nanochannel shown in Fig. Ib. In this case the total 
system energy E is given by a sum over all Coulomb 
interactions. Due to strong nonlinearity of the system 
the minimal energy configurations should be find numer- 
ically using the methods described in [17, 20-22]. We 
take a finite number of electrons N for L periods of a 
channel of finite length. In numerical simulations we put 
the channel on a cylindrical surface in 3D with electron 
coordinates being x\ = Lsm(si/L), j/j = L cos(si/ L), 
z = asin(sj) where Si is coordinate along channel for 
electron i. Thus the channel, filled by N electrons, wig- 
gles in z— direction making L periodic oscillations along 
cylinder of radius L with periodic boundary conditions. 
The Coulomb energy of the system is 



£ = ]ri/i?( Sl , Sj ) 



(i) 



j>i 



where R(si, Sj) is the distance between two electrons. We 
find from geometry R 2 (si,Sj) — 4L 2 sin 2 [(sj - Sj)/2L] + 
a 2 (sin Si — sinsj) 2 . Here we choose dimcnsionlcss units 
for charge e and length, so that the channel period length 
is £ = 2tt and dimensionless amplitude of channel oscil- 
lations is a. In the limit of large L we have a channel 
wiggling in (x, z) plane. At a — we have electrons on a 
circle. 

The equilibrium static configurations are defined by 
the condition dE/dsi = with a minimal ground state 
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FIG. 3: (Color online) Dependence of the dimensionless 
phonon gap A/e c on the electron density v — N/L for 
a = 0.7(a), 1(6), 1.2(c), 1.5(d). Here L = 89 (black), 233 
(gray/red). The straight line shows empiricai dependence 
A/e c oc (N/L) 1/2 for (c,d), where e c = 2ne 2 v/t = v is the 
Coulomb energy. 



energy configuration determined numerically by the stan- 
dard methods [17, 20, 21]. Using these methods we also 
find the phonon Spectrum co(k) of small oscillations at 
the ground state. It is easy to see that the total en- 
ergy E is invariant for a homogeneous shift of all elec- 
trons by Ss when the distance between nearby electrons 
is Sj-|_i — Si = 2itm that corresponds to electron den- 
sity v = N/L with resonant rational values v m = l/m. 
Hence, at such a density the Wigner crystal can freely 
slide along the channel. For irrational density values the 
properties of sliding are much more tricky. An example is 
shown in Fig. 2 for two very close incommensurate den- 
sities v. For v = 239/233 we have a smooth hull function 
Si = h(xi)(mod2Tr) with the gapless phonon Spectrum 
lo oc ujQk/N at small wave numbers k/N. Here Xi(mod2ir) 
are ground state electron positions at a = 0. The di- 
mensional unit of frequency luq = (e 2 /(m(^/27r) 3 )) 1 / 2 is 
expressed via the particle charge e, mass m and channel 
period i, we omit it in the further dimensionless com- 
putations. This regime corresponds to the continues in- 
variant Kolmogorov-Arnold-Moser (KAM) curves as it is 
discussed for the Frenkel-Kontorova model [17-21]. Here, 
the Wigner crystal can slide freely along the nanochan- 
nel. In contrast, for very close density v = 244/233 the 
hull function starts to have devil's staircase form, well 
known for the cantori regime in the Frenkel-Kontorova 
model. Here, the gap A appears in the phonon Spectrum 
so that the crystal is pinned in the channel. This regime 
corresponds to the insulating phase. 

The dependence of phonon gap A on electron density 
v is shown in Fig. 3 for various values of channel defor- 
mation amplitude a. At small deformations the gap is 
zero for a large fraction of densities v (Fig. 3a) and the 
crystal can slide freely along the channel. However, at 
larger deformations the gap disappears only in a vicinity 
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FIG. 4: (Color online) Dependence of rescaled phonon gap 
A/e c on channel deformation amplitude a at various values 
of electron density v with the number of electrons N = 241 
(black), 269(blue), 337 (magenta), 377 (red), 307 (green) 
(curves from right to left at A/e c = 0.01 respectively) at 
L = 233. 



of rational densities v m (Fig. 3b, c) and at strong deforma- 
tion regime only narrow zero gap intervals remain around 
these density values (Fig. 3d). We note that our numer- 
ical data are obtained at rather large number of elec- 
trons N and channel periods L so that the dependence 
A(i/) found numerically corresponds to the limit of infi- 
nite channel length. Indeed, the function A(i/) remains 
practically unchanged with an increase of L (Fig. 3c). 
The global dependence of A on v corresponds to fre- 
quency of small charge oscillations A oc v z / 2 oc l/£ 3 ^ 2 , 
being in agreement with data of Fig. 3c, d. 

A remarkable feature of the dependence A(i/) is its 
very sharp variation with density v and deformation o. 
The dependence is enormously sharp in a vicinity v = 1: 
for v < 1 there is crystal sliding in the channel while 
only slightly above v = 1, e.g. for N/L = 1 + 11/233, 
we obtain the insulating phase. This reminds a sharp 
changc of conduction propertics of organic materials with 
pressure [11] which effcctivcly modifies v and a values. 

The dependence of phonon gap A on channel defor- 
mation a is shown in Fig. 4 for few density values v. 
The gap changes smoothly with a for a > a c where a c is 
a certain critical value of deformation which depends of 
density v. For a < a c we find very sharp drop of A which 
becomes exponentially small, e.g. A drops by 5 orders 
of magnitude when a is decreased by a couple percents 
in a vicinity of a c . Since simulations are done at a finite 
N this means that in the thermodynamic limit A = 
for a < a c . We interpret these data in a way similar 
to the case of the Frenkel-Kontorova model [18, 20, 21]: 
for a < a c (u) we have an analytic invariant KAM curve 
with a rotation number corresponding to a given density, 
while for a > a c this curve is replaced by a cantori with 
a finite phonon gap and pinning of the crystal. 

To understand in a better way the numerical results 
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FIG. 5: (Color online) Map kick functions g^((j>) (a) and g v (v) 
(b) obtained from the groundstate electron positions s, in 
nanochannel (points), full red/gray curve in (a) shows the 
theoretical dependence from the map (2), see also text. Here 
N=377, L=233, a = 0.5. 

presented above we derive an approximate dynamical 
map which determines recursively the electron positions 
along the channel. The recursion is given by equilibrium 
conditions dE/dsi = 0. Assuming that o<C lwe can ex- 
pand R in a that, after keeping only nearest electron in- 
teractions, gives recursive relations between Sj_i, Sj, Sj+i. 
They can be presented in a form of dynamical map 



v + 2a 2 (1 — cos v) sin 2<f> 
4> + v + a 2 sin v cos 2<j> , 



(2) 



where v = Sj — Si—i, <f> — s, are conjugated action-phasc 
variables, bar marks their values after iteration. The map 
is implicit but symplectic (see e.g. [23]). To check its va- 
lidity we use the values s, obtained for the groundstate 
configuration and extract from them the kick function 
gø = sin 2ø from the values v — v — 2a 2 g v {v)g ( j 1 {4>) with 
g v (v) = 1— cosi;. Such a check shows that the map (2) in- 
deed gives a good description of actual electron positions 
Si up to moderate values of a (Fig. 5). 

At small a the phase space of the map is covered by 
invariant KAM curves as it is shown in Fig. 6 (left). At 
larger a a single chaotic component covers a significant 
part of the phase space (Fig. 6 right). Locally the dy- 
namics is approximately described by the Chirikov stan- 
dard map with the chaos parameter K ps 4a 2 (1 — cosv). 
According to [19, 23] the KAM curves are destroyed at 
K > 1 that is in a good agreement with our numerical 
data of Fig. 3 where the KAM curve with the golden 
rotation number v — 0.618... goes to the cantori regime 
approximately at a w 0.4. We note that at small charge 
density v the parameter K is small K » 2c? v 2 <§C 1 that 
corresponds to the KAM regime and a conducting phase 
of Wigner crystal in agreement with the data of Fig. 3. 

Of course, the map description is valid only up to mod- 
erate a values. At a > 1 the expansion in a is no more 
valid that explains the asymmetry in the dependence for 
A(v) at v < 1 and v > 1 which is absent in the ap- 
proximate map (2) but is clearly seen in Fig. 3. Furthcr 
studies are required to obtain a map description at large 
values of deformation a. 
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FIG. 6: (Color online) Poincare section for the dynamical 
map (2) at a = 0.25 (left panel), 0.5 (right panel). 



Our studies determined conditions of sliding and pin- 
ning of the Wigner crystal in snaked nanochannels. Hcrc, 
we performed only classical analysis. According to results 
of [17] the quantum effects are weak if the dimension- 
less effective Planck constant % e ff — (27r7z 2 /me 2 ^) 1 ^ 2 is 
small. This is defmitely the case for supercapacitors with 
£/2tt ~ lnm, large ion mass m ~ 4 • 10 4 m e compared to 
electron mass m e , that gives h e ff ~ 10~ 3 . The charge 
storage process in this case starts with small charge den- 
sity values v where the ions slide easily in nanochannels 
since the gap A is practically absent there (see Fig. 3). 
However, with the increase of v ions form the Wigner 
crystal which is pinned inside the nanopores at large 
v values. We think that this is the physical mecha- 
nism behind the charge process of electrochemical ca- 
pacitors studied in [13, 14]. We note that the energy 
of pinned Wigner crystal can be estimated as Ww ~ 
Sde 2 /e(£/2n) 4: , where S is the surface area, d is the deep- 
ness of nanopore layer on the surface and e is the dielec- 
tric constant. For typical parameters e = 5, £/2n = lnm, 
d = 1/iTO we obtain Ww /S ~ 5- 10~ 3 J/cm 2 . It is natural 
to assume that a part of this energy can be used during 
recharging process that makes it comparable with the 
surface energy density reached in supercapacitors with 
W/S ~ 10~ 3 J/cm 2 at maximal capacitance per area 
C fa 400^F/cm 2 and voltage U ~ 2V [13, 14]. We note 
that our estimate gives an increase of Ww with a de- 
crease of nanopore size i that qualitatively corresponds 
to the behavior observed experimentally (see e.g. Fig. 3a 
in [14]). At the above parameters the typical pinning frc- 
quency is woA/27r ~ 50GHz so that the Wigner crystal 
should be sensitive to microwave radiation in this fre- 
quency range. 

In contrast, for CDW in organic conductors [11] we 
have m ~ m e , £/2ir ~ 3A that gives h e ff ~ 0.5 so 
that quantum effects can play an important role. Fur- 
ther studies are required to analysis quantum properties 
of crystal sliding but we expect that they will have sim- 
ilarities with the quantum Wigner crystal in a periodic 
potential [17] and the quantum Frenkel-Kontorova model 
[22] . The classical pinned phase should correspond to the 
insulator phase, while we expect that the classical sliding 



phase may correspond to the superconducting regime in 
the quantum case. Indeed, the sliding phase has a linear 
dispersion law w(k) that can be similar to the situation in 
superfluid phase. The sharp transitions from conducting 
to insulating phase with charge density variation are well 
pronounced in the classical regime and are expected to 
be present also in the quantum case. This can be at the 
origin of high sensitivity of organic conductors to pres- 
sure. Further studies should shed more light on the quan- 
tum properties of Wigner crystal in snaked nanochannels 
and organic molecular chains. It would be very interest- 
ing to study such effects experimentally creating artificial 
snaked channels with electrons on a surface of liquid he- 
lium [9]. 

This research is supported in part by ANR PNANO 
project NANOTERRA. 
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